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Abstract 



> 

o, 

^^ , We propose a general algebraic analytic scheme for the spectral trans- 

form of solutions of nonlinear evolution equations. This allows us to give 

>— P ' the general integrable evolution corresponding to an arbitrary time and 

space dependence of the spectral transform (in general nonlinear and with 
non-analytic dispersion relations) . The main theorem is that the compat- 
iblity conditions gives always a true nonlinear evolution because it can al- 
ways be written as an identity between polynomials in the spectral variable 
k. This general result is then used to obtain first a method to generate a 
new class of solutions to the nonlinear Schrodinger equation, and second to 
O ' construct the spectral transform theory for solving initial-boundary value 

^ ' problems for resonant wave-coupling processes (like self-induced trans- 

parency in two-level media, or stimulated Brillouin scattering of plasma 

, y waves or else stimulated Raman scattering in nonlinear optics etc.). 
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1 Introduction 

The spectral transform method to build and solve classes of nonlinear evolution 
equations (for a field Q{x,t)) works as a nonlinear extension of the Fourier 
transform |^] : it associates to the field Q{x, t) its spectral transform R{k, x, t) 
which is a distribution in the spectral complex variable k and a function of the 
real space and time variables x and t. The domain of definition for k and the 
{x, i)-dependence of R uniquely determine both the related nonlinear evolution 
for Q and the class of functions to which it belongs. 

For instance when k varies on C and when i? is a 2 x 2 off-diagonal matrix 
depending on x and t through 

i?, = [R, ifcfTg], Rt = [R. ik'cTs]. (1.1) 

then Q is a 2 X 2 off-diagonal matrix obeying the nonlinear Schrodinger equation 

«^3Qt = -^Q... + Q'- (1-2) 

This is the very reason why the spectral transform allows us to solve the nonlin- 



ear evolution (1.2) : the corresponding evolution in the spectral space is linear. 
It is then enough to have a bijection between Q and i?, which is the set of direct 
and inverse spectral problems. 

The main purpose of this work is to answer the following simple question: 
can we still build and solve nonlinear evolutions when the (x, i)-dependence 
of the spectral transform is arbitrary. We shall prove in particular that the 
nonlinear evolutions associated with a local 2x2 matrix d problem on the 
Riemann sphere can always be written in closed form whatever be the space 
and time dependences of R in a very general class of equations (not necessarily 
linear). 

Here the d problem will be understood as being the main tool which links 
Q{x, t) to R{k, X, t) by means of a matrix valued function /x(A:, x, t) solution of 

M(fc) = 1 + — yy ___^(A)i?(A) (1.3) 

(being parametrically induced by R{k, x, t), the (x, t)-dependences can be omit- 
ted everywhere in /x and R). If we operate with the 9-operator (see next section) 
on both side of the above equation we obtain the so-called d problem for ^ : 

t-^^ik) = fi{k)Rik), Mfc),^ 1 (1.4) 

ok fc^oo 

It is clear that any dependence of R(k) on external independent (real) parame- 
ters {x and t) induces a dependence of /i(fc) generally through an overdetermined 



system of differential equations : the so-called Lax pair. Then this overdetermi- 
nation implies some constraints on the coefficients of the Laurent series for /Lt(fc) 
(see ([2.15) next section). When they form a closed system independent of the 



variable k, these constraints are nothing but the nonlinear evolution equation 
related to the given space and time dependence of R{k). 

The result presented here will undercover some interesting properties of 
overdetermined s yste ms (for R(k) and for /i(fc)) when they are related through 
a B problem like (|l.4|) . First we shall discover that the differential system for R 
(data of Rt and Rx) although very general, has to be written in a special form 
in order to obtain a non-trivial integrable related evolution. This remark will 
actually bring into light the necessity of a spectral parameter in the Lax pair. 
Second, our result will be shown to generate new classes of integrable systems 
with in particular some very interesting applications to problems of interaction 
of radiation with matter (for instance stimulated Brillouin and Raman scat- 
tering), which are always associated with a mixed initial and boundary value 
problem. 

Although some subcases and applications of the material presented here have 
been reported elsewhere in P|...|q], the results presented here are mainly new 
and we shall avoid using directly any previous result. In that purpose, the next 
section is a brief recall of the basic mathematical tools that we have chosen to 
present in their light version, reporting to M for mathematical details. 

The section 3 is the heart of this paper : the main theorem is established 
and we obtain the general algebraic-analytic scheme of the spectral transform 
method. Some important consequences of general interest are then discussed. 

In section 4 we present an application of the general theory to the nonlinear 
Schrodinger equation (NLS) and exhibit in particular a new method to generate 
a more general class of solutions. 

The section 5 is devoted to the main interesting application of the general 
result: the theory is used to build integrable coupled systems for which some 
arbitrary boundary values are prescribed. As a consequence the time evolution 
of the spectral transform is in general nonlinear with some very new properties 
of the solution in physically interesting cases. 



2 Notations, definitions and basic tools 

The derivative with respect to fc,(/c = kn + ikj) formaUy defined as 

d__l(_d_ ._d_\ 
dk~ 2\dkR^'dki) ^ ' 

vanishes on the space of analytic functions of k (it is sometimes cahed the 
measure of non-analyticity). It can be conveniently defined by the Cauchy- 
Green representation of a function in C: 

/(fc) = /o(fc) + /i(fc) 

sm = ^jf^Ji^d_m (2.2) 

■'^^ ' 2nrJJ X-k dX ^ ' 

where /o(fc) is the asymptotic behavior (analytic in fc) of /(fc) when k -^ oo. 

Here we shall work more particularily in the space of function f{k) whose 
behavior fo{k) is a polynomial (in k) which we note Vk, that is 

Mk) = Vkif} (2.3) 



We call Tip the space of functions which have the representation (2.2) with the 



constraint (2.3) 



A consequence of (2.2) is that we can formally write 

d 1 



= 2iTTS{X - k) (2.4) 



dkX-k 
where 5 is the Dirac distribution in C, namely 

dX A dXf{X)S{X -k) = f{k) 
dXR dXif{Xn + iXi)-S{XR~kR)5{Xi-kj). (2.5) 

-oo J— oo ^ 

We shall also use the distributions S"^ defined by 

/+00 
dXRfiXR±iO), (2.6) 

oc 

and make use of the notation 

f^{kR)^f{kR±tQ). (2.7) 

From ( ^.4[ ) we can deduce 



d f+°° dXr 



., , « -f{XR) = -7:f{k)d^ikj). (2.8) 

(jk j_oo Xr — k ± tu 



Let us recall also the Sokhotski-Plemelj formula 



1\ p-\-'DC 1\ 

-f{X) = ±^7rf{k) + P ^^/(A), (2.9) 



, \-k±iO-'^ ' ■' ' ' J_^ X~k 

whereijf stands for the Cauchy principal value integral. The two above formulae 
allow us to obtain by direct calculation 

'- [f{k)S+{kj) .f{k)S-{kj)] = tf{k), (2.10) 

which provides the d derivative of a function discontinuous on the real axis but 
analytic elsewhere. 

We note that the Cauchy-Green integral equation ( ^.2| ) can also be used to 
define the inverse of the d operator in Tip as 

d-'{^^)=nk)-v,{f}, (2.11) 

which is nothing but fi{k). Hence the k-polynomial part of /(fc) appears as the 
constant of integration in Tip of the d operator. Whenever possible we use the 
short-hand notation 

«/^|/'^>' S-/^i//^/(A) (2.12, 

and we have the very important property 

a-i/ = 0(i) as fc->oo (2.13) 

Next we adopt for the Laurent series of f{k) the followmg notation 

oc 

/(fc) = /o(fc) + ^fc-"/<"\ (2.14) 



where from (2.2) 



Jin) ^Zl ffdXAdX ^B. A"-i (2.15) 

2t7TjJ dX ^ ' 

For a complex valued function /(fc) = a{k) + ib{k) of the complex variable 
k = kji + iki, we note 

/(fc) = a{kR + iki) - ib{kii + iki), (2.16) 

/*(fc) = aikji - iki) - ibikji - iki) = J{k). (2.17) 



Then, from the definition ( |2.6D , since S^ are real valued distributions obeying 
S~^{ki) = 6^{—ki), we have 

6+{kj) = {6-{kj)r. (2.18) 

It is also usefuU to remember that, when A depends on a real parameter, say 
t, we may write 

d c .^ 5(k — A) d\ 

which is meaningfuU on the space of functions of k that vanish in k ~ X. 

Finally, we shall work throughout the paper in the space of complex-valued 
2x2 matrices, say M, which elements rriij are differentiable functions of two real 
variables (x and t) and distributions in a complex variable (k). It is convenient 
to split the matrices in their diagonal and anti-diagonal parts, namely 

rD , n,rA nrD _ I "^11 \ n^A f "112 



M==M^ + M^, M"" ^ [ '' , M-^^ n ]■ (2.20) 

V TO22/ \m2i J 

Using the standard definition of the Pauli matrices 

we have for instance 

[cr3,M] = 2o-3M^. (2.22) 



3 General theorem 



The starting po int is the Cauchy-Green integral equation (1.3) or equivalently 



the 9-problein (L4) for the matrix valued complex function /i(fc) : 



^^{k) = ^iik)Rik), Mfc) = l + 0(i), (3.1) 

where the given datum R{k) is a distribution in the set of 2 x 2 matrices of 
vanishing trace: 

tr{R{k)} = 0. (3.2) 

Theorem 1: 



When R{k) depends on two real parameters x and t according to 




i?:, = [R, A] + N, 


(3.3) 


Rt = [R, ft] + M, 


(3.4) 


the quantities 




Uo = -n-i/iA^-i}, Ui = d-'{fi{N - dA)^i-^} 


(3.5) 


Vo = -Pkifinfi-'}, Vi = d-'{fi{M - dn)^i-'} 


(3.6) 


obey the following equation 




Uo,t - Vo,^ + [Uo, Vo] + VkifiAfi-'} = Vk{[Vo, U,] - [Uo, V]}, 


(3.7) 


where the matrix A stands for 




A:^ At-n^ + [n,A], 


(3.8) 


and where the entries {A, N, ft, AI} obey the compatibility constrain 


b 


Nt - [N, 0] - M^ + [M, A] = [A, R]. 


(3.9) 



Equation (3.7) is a nonlinear evolution equation because it is an identity be- 
tween polynomials in k and then it has to be read as a set of local fc-independent 
differential equations (each coefhcient of a given power of k vanishes) . 

Here above the 2x2 matrices A and fl are functions in Tip (see sec. 2) and M 
and N are in general distributions. All these 4 given matrices are differentiable 



fuctions of the real variables x and t. Up to now the {x, i)-dependence (3.3)(3.4) 
of R is quite arbitrary but we will understand later why this specific structure is 
important (rather than just letting A := fJ = 0). Note also that nothing prevents 
these equations to be nonlinear in R. 



The problem is to find the {x, i)-dependence of the solution fi{k) of ( |3.1[) 
when R{k) is required to solve ( |3.3[ ) and (3.4). The method proposed in W 
consists in using first the commutativity of the differential operators. 



d_ d_ 
dt' dk 



= 0, 



d_ d_ 
dx' dk 



= 0, 



(3.10) 



and second the property that /^(fc) is a regular matrix. Indeed from ( [3.2| ) we 
easily obtain 

d 

-^det^ik) = (3.11) 

ok 

and hence, since fi{k) = 1 + 0(-t:) for large k, the Liouville theorem implies 



detfi{k) 



(3.12) 



These results allow us to obtain by direct computation the following relations 
d{ifi,-iiA)ii-^}^fi{N-dA)^i'^ (3.13) 

d{{iit - M^)At"^} = KM - dn)ii-^ (3.14) 



They can be integrated by using (2.11), the result can be written as the following 
overdetermined system 



fix = ^J.A + U^l , fit = /ifi + Vfi, 



in which U and V are in Tip, that is 



U = Uq + Ui , y^T/^o + V^i, 



(3.15) 



(3.16) 



where Uo, Ui, Vo, Vi are given in (3.5) and (3.6). 



Consequently the d problem (3.1) serves as the link between the two sets of 
overdetermined differential systems (3.3) (3.4) and (3.15). Computing R^t in 
two ways leads to the constraint ( |3.9D and doing the same for ii^t leads to the 
compatibility condition: 



Ut-Vt + [U, V] = -fiAfi- 



(3.17) 



We come now to the demonstrat ion o f (3.7). This is done fi rst b y co mpu ting 
the quantity Ut — Vx appearing in ( 3.17 ) with the definitions (|3.5|) to (3.6) . To 
arrange conveniently this quantity we have to use the inverse of (3.5) and (3.6) 
that is 

dU ^fi{N -dA)fi-^ , BV ^ ii{M - dn)^-^ (3.18) 

with this help one can write 

Ui,t-Vi,x = d-'{fiiNt-Mx-[N,n] + [M,A])^-'~ 

- fid{At-nx + [n,A])fr^-[u,dv]^[du,v]}. (3.19) 



At this stage enters the constramt (3.9) which guarantees the compatibility 
of the choices for Rx and Rt ■ Then the above equation can be written 



Ui,t - Vi,. 



-d-^MdA + [R, A])^i'^ + d[U, V]} 



or else 



Ui,t - V^i,. = -d-'d{f,Af,-' + [U, V]} 



(3.20) 
(3.21) 



The inversion formula ( 2.11 ) then gives 

Ui,t - Vi,x = -fiAfi-^ - [U, V] + Vkit^A^i-' + [U, V]} (3.22) 

which then, inserted in (3.17), easily leads to ( |3.7D fo we have obviously 

Pk{[U, V]} = [Uo, Vo] + Vk{[Uo, Vi] + [Ui, Va]}. (3.23) 

First we note that this result is far from being trivial. Indeed, looking at 
the structure of U and V given in (|3.5D and (3.6), we realize that the equation 



(3.17) involves quite complicated functions of k through many d ^ operators 



(see ( ^.12 )). The relative simplicity of the result (3.7) originates from the fact 
that all terms in the integral part of (3.17) can actually be written as an exact 
9-derivative. 

Second it is the very reason why the spectral transform works also in the 
case of singular despersion relation. In previous works (f.i. in H to H) the k- 
dependence in ( |3.17 ) was seen to cancel out as a result of the particular choice 
of A but in no way as a very general property. 



Third the result ( |3.7| ) justifies the choice of the structure in ( |3.3[ ) (3.4). 
Indeed, avoiding any structure in (3^) and (3^) by just lettin g il = A — , 
makes the equation (3/7) tri vially solved: we would have from ( |3.5| ) and (3.6) 
Uq = 0,Vo = and hence ( 3.17 ) would not lead to any nonlinear evolution 
equation. 



Note also that another necessary condition that ( 3.17 ) be a n onlinear evo- 
lution equation is that Uq be nonzero, that is from (|3.5|) and (3.1) that ■PfcjA} 
be non zero. This property can be un derstood as the necessity of a spectral 
parameter in the spectral problem ( 3.15 ) (here PfejA} would play the role of the 
spectral parameter). 

We have finally built a set of relations between partial differential equations 
(obtained by means of the 9-problem) which has to be understood as follows: 
if the set {i?, $7, A, M, iV} obeys the system (b3[)(3.4) and the compatibility 
constraint (|3.9[), then the set {f/o, C/i, Vb, Vi}, given from the solution fi of the 



9-problem in (3.5) (3.6), obeys the nonlinear evolution (3.7) 



The compatibility constraint ( p.9[ ) will be seen to be essential. At this stage it 
is important to note that the degree of complexity of this constraint is usually 



greater than that of the nonlinear equation itself. Hence the main point to 
make here is that, in order to build from theorem 1 some practical tools, we 
have to make a decision, for instance by choosing an elementary solution of the 
constraint (|3.9|) 
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4 The nonlinear Schrodinger equation 

A n atural choice for the arbitrary functions A and Q, for which the constraint 
( p.9| ) is independent of the spectral transform R, is 

At-n^ + [fl,A]=0. (4.1) 

It is usefuU to rewrite the theorem in this subcase: 
Theorem 2 



For 




R^ = [R, A] +N, Rt = [i?, n] + M, 


(4.2) 


At-n, + [n,A] = o, 


(4.3) 


we have the nonlinear evolution equation 




Uo,t - Vo,^ + [C/o, Vo] = Vk{[Vo, Ui] - [Uo, Vi]}, 


(4.4) 


with 




Uo = -T'fciMM"'}, c/i = a-i{M(iv - ^A)^l~'}, 




Vo = -T'fel/if^M"^}, 14 = a-i{M(M - ar!)/i-i}, 




if M and A^ satisfy the constraint 




Nt-[N,n]-M.,. + [M,A]=0. 


(4.5) 



The simplest non-trivial choice is realized by taking for A and Q two diagonal 
matrices, polynomial in k with constant coefficients. In that case, the matrix A 
given in ( |3.8D does vanish and a trivial solution to (p^) is Af = A^ = 0. Then 
we recover the usual hierarchies of integrable partial differential equations. To 
set an example, let us chose 

A^ikas, n^ik^as, M^N = 0. (4.6) 

Then from (jsj) and (|U) 

Ui =0, Vi^ 0, (4.7) 



and from the Laurent series (2.14) for /i(fc) 

C/o--^fca3 + ^[a3,/i^'^] (4.8) 

Vo = ^lasP + ik[a3, Ai''^] - ^[a3, /i^' V^'^ + «[^3, M^'^] (4-9) 

The above two polynomials in k are then inserted in the evolution equation 



(4.4) and each coefficient of any power of k is required to vanish. Defining the 



anti-diagonal matrix Q{x,t) by 

Q = t[a3,ti^% (4.10) 
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we easily obtain: at order fc^ the equation (4.4) is automatically verified, at 
order k this equation gives 



a3,g/i^i^-*[a3,Ai^2)] 



^X • 



(4.11) 



The above equation can the be shown to be equivalent to the diagonal part of 
the fc-independent term in (4.4), and it actually defines the quantity i[(T3,/x'-2)] 
from Q. Provided this information, the anti-diagonal part of (4.4) finally gives, 
again by using (4.11), the nonlinear Schrodingcr equation 



la^L 



0. 



(4.12) 



At this point there is the important remark that A and Vt can be modified in 
such a way as to keep the evolution equation unchanged (that is to keep C/q ^^nd 
Vo unchanged) but to change the (x, t)-dependence of R{k,x,t). This will then 
generate, through the solution of ( |1.4| ) and of (4.10), a larger class of solutions 
of the NLS equation. This can be done by choosing 



Vk{A} = ikas, Vki^} = i/cVg, 



(4.13) 



N = dA, M = dfl, 



for which Uq is still given by 



(4.14) 
and Vq by (O), and for which Ui ~ Vi ^ 0. 



Note that the constraint ( |4.5| ) is automatically verified from (4^) and ( 4.14 ). 
While R{k, x, t) now obeys a much more general (x, i)-dependence, the re- 



lated nonlinear evolution is still the above NLS equation (4.12), and hence we 
have indeed here a method to generate a larger class of solutions to the NLS 
equation. The method to perform this can be sketched as follows: 
Step 1: get some explicit functions A and fi by means of 



At-nx + [n,A] = 

Pk{A} = ika3, VkM = ik'^as 



{A{k,x,t),n{k,x,t)}. 



Step 2: solve the (a;,i)-dependences for R{k,x,t) 

-^{R{k,x,t)}. 



Rx 
Rt 



[R,A] 



dA 

on 



(4.15) 



(4.16) 



Step 3: solve the linear Cauchy integral equation for fj,{k,x,t) 

[l^ik)^l + i^Jjq^fi{\)R{\)]^Mk,x,t)}. (4.17) 

Step 4: compute the solution Q{x, t) of the NLS equation by means of fi{k, x, t) 

Q{x,t)^i[a3,^i'^'Hx,t)]. (4.18) 
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Steps 3 and 4 hereabove are the usual problems to be solved when a solution 
to an integrable nonlinear evolution is seeked and step 2 is a linear differential 
equation which is solved through standard methods. The novelty here lies in the 
step 1 which solution actually reduces to solving the nonlinear integrable hier- 
archy generated through the potentials A and 51 by the compatibility condition 
( [l^ ) . The problem is that of the stringent constraints 

A^ika-i + Oih, n = ik^a3 + Oil;). (4.19) 

The technical complexity of the proposed method is such that we have not 
been able to find an explicit example. Note in particular that, when A is not 
simply ika^, the scattering problem (the x-part of the Lax pair (3.15|)) has to 



be rebuilt completely. And this is a necessary step in order to find the constants 
of integration for the two differential equations in (4.16|). 



In this context, the last point to check is the unicity of the solution of 
the Cauchy problem for the NLS eq. when Q{x,0) is in the Schwartz space. 
This unicity will hold if, when analyzing Q{x, 0) through the Zakharov-Shabat 
spectral problem, in other words with 

A == ikas, (4.20) 

we would obtain the usual time evolution of the spectral transform, in other 
words 

n = ik^as. (4.21) 

This result is obtained here by first remarking that the compatibility equation 
( [4.3| ) in which A = ika^ gives by induction (successive powers of 1/fc) 

fl"^ ^ 0. (4.22) 

Second, we use the property given in the appendix that the spectral transform 



R is an anti-diagonal matrix. Then the time evolution in (4.2) can be splitted 
in its diagonal and anti-diagonal parts to give 



dn"" = [n"", i?], Rt = [R, n^] + an"". (4.23) 



The first of the above equations can be explicitely integrated thanks to (4.19) 



and (4.22) and we have 

n^ = ifcVg (4.24) 



which acheives the proof. Indeed, (4.23) now reads 

Rt = [R,ik''a3]. (4.25) 
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5 Integrable asymptotic boundary value prob- 
lems 

5.1 Statement of the problem 

The general theorem 1 is shown here to serve as the basic tool for demonstrating 
the following important result. 

Theorem 3: 



The syste 


m of 


coupled equations for the three fields 


q{x,t), fli 


(k, X, t) and 


a2{k,x,t) 




qt = 

a2.,x 


r+co 

/ gdk ai02, 

ai,x = qa2, 
— 2ika2 — fjqai 




(5.1) 


(with a = 


= ±, 


x,k G R, t 


> and g = g{k) an 


arbitrary 


function in 


L"^ which 


may 


also depend on t), is integrable for the arbitrary asymptotic 


boundary 


values Oj(fc, ±oo,i) 









More precisely we will give the evolution of the spectral transform in the 
following four cases: 

ai -^ h{k,t), 02 -^ /2(fc,i)e2^'=^ (5.2) 



2ikx 

a:— i- + cxD a;^ — oo 

2ikx 

a;— > — oo a;^ — oo 

ai ^ Li{k,t), 02 -^ L2{k,t)e''''\ (5.5) 



ai -^ Ji(fc,t), 02 -^ J2{k,t)e^'''^ (5.3) 

a:— i- + oo a;^ — oo 

ai -^ Ki{k,t), a2 -^ K2{k,t)e^'''^ (5.4) 



Here above the asymptotic boundary values {Ii, I2, Ji---) are arbitrary functions 
of the real variable k in L^ with an arbitrary dependence on the real parameter 



t. Of course, due to the linearity of the equation for Oj in (5.1 ), any of the above 
behaviors can actually be obtained from a particular one. However we will see 
that it is quite usefuU to have the four evolutions of the spectral transform 
separately as the link between these cases can well be very complicated. Note 
also that to each of the above 4 cases there are two choices a — + and a = — 
for which the corresponding evolutions of the spectral transform will be seen to 
be quite different. 
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The above system will be understood as the equation for the two envelopes 
ai and 02 of the two components of a high frequency field, interacting resonantly 
with a low frequency field of envelope q. The parameter k actually measures a 
frequency missmatch and g{k) measures the relative intensity of the coupling of 
the waves for each value of k. Hence g{k) has to be positive, symmetric around 
the proper resonance fc = (even function of k) and of finite support. Conse- 
quently the system (|5.l| ) is in three variables and is constituted of a differential 
eq. in the variable x running on the infinite line, a differential eq. in the variable 
t running on an arbitrary compact support and an integral eq. in the variable k 
running on a given compact support (the support of the measure g{k)). Given 
g{k), g(a;,0) and the boundary values for the a^-'s, the system is closed and the 
problem of solving it is well posed. 

The proof of the integrability of ( 5J_ ) has been first given in [^ Q and follows 



a large series of papers devoted to the self-induced transparency (SIT) equations 
of McCall and Hahn |^ which were given a Lax pair and soliton solutions by 
Lamb P|. The sharp line case was then shown to be completely integrable in 
[ [lO| where the N-soliton solution is displayed. Then the inhomogeneously broad- 
ened case was also shown to be completely integrable in [[111 where in particular 
the evolution of the continuous part of the spectrum has been obtained, which 
has allowed the authors to demonstrate mathematically that this simple model 
does contain the property of transparency (the radiation was shown to vanish 
exponentially as the input laser pulse propagates in the medium). Later, more 



general asymptotic boundary conditions (a subcase of (5.2) for |/i|^ -I- I/2P = 1) 
have been studied in |l3] and shown to lead to interesting physical consequences. 
Simultaneously the scientific community became aware of the importance of 
this system as a fundamental model for resonant wave coupling processes in the 
slowly varying amplitude approximation (SVEA), role which is played by the 
nonlinear S chro dinger equation in the scalar case (no wave coupling). Indeed, 



the system (5A) appears as a universal limit in many different physical systems, 
not only for the Maxwell-Bloch system in two-level one-dimensional media (the 
starting equations for SIT). It appears for instance in nonlinear optics to de- 
scribe stimulated Raman scattering [ [l3[ , in diatomic chains of coupled harmonic 
oscillators [Q, in plasma physics for stimulated Brillouin scattering ||l^... The 
point is that to each individual physical situation, there correspond different 



asymptotic boundary values ||6|. Consequently, the integrability of (5.1) has 
not a unique common implication on the behaviour of the solution and, as we 
shall see, the properties of the solution can well be drastically different from a 
problem to another. 



The integrability of the system (5.1) is proved in three steps. First we obtain 



a general integrable system of a structure similar to (5.1) for a particular choice 
of the entries (TV, A, M, $7) in theorem 1 (sec. 3). Second we consider some 
convenient reductions of the obtained evolution equation. Third the values of 
the arbitrary functions occuring in the chosen set (iV, A, M, i7) are determined 
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in terms of the asymptotic boundary values. The proof wiU actually only be 
sketched for the details can be found in IIgII. 



5.2 General integrable system 

Let us choose 



M 



N = 0, A = ik(j3 
m-{k,t) 



eKp[2ika3x] 



m+{k,t) 

dn=p{k,t)a3, ^-{^1 = 



(5.6) 

(5.7) 
(5.8) 



For that simple choice the compatibility constraint (3.9) ho lds (note that A — 0) 
and the potentials U and V can be computed out of ( |3.5| )(3.6) and read 



f/o = -^fca3 + ^[a3,M'l'], C^i = 0, (5.9) 

^0 = 0, Vi=d-^{^i{M ^dVCj^i-^}. (5.10) 

The 2x2 matrix /i(/c, x, t) obeys the differential equation ( |3.15|) namely here: 



-ifc[cr3,^]+ ( ^ ^ ) /i, 



where we have defined 



q 
r 



^[cr3,^l 



1)1 



(5.11) 
(5.12) 



Finally the evolution equation (3.7) of theorem 1 can be written 
1 



qt = 



n 



dk A dfc(2^ii^i2 p^- ^ii^m e""^'"^ - /^a "^ e^'"^), 



(5.13) 



dk A dk{2fi2if^22 P + A^2i ^^ ^ 



-2ikx 



/4 ^^ e'-'n, 



where of course the /i^j 's are the matrix elements of /i(fc, x, t). 

The equations (5.11) and ( ^.13| ) constitute a system of 6 coupled fields to 
which we have to associate the initial data (^(x, 0), r{x, 0)) and some boundary 
values for the /Xy's (for instance aX x = cx)) for all t. 
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5.3 Constraints 



We shall work out three stages of simplifications of the system (5.11) (5.13). 
The first one is a result of the constraint on the space of funtions to which 
belongs the initial datum (we call it a structural constraint) . The second one is 
a constraint on the nature of the coupling term in the evolution equation to be 
solved (we call it a natural constraint) . The last one is a reduction of the 6-field 
system to a 3-field one by assuming a simple relation between q and r (we call 
it a reduction constraint). 

i) Structural constraints 



To ensure consistency of the evolution (3.4) as an equation for distributions, we 



must impose the constraints which result from the support of the distribution 
R{k, X, t) in the complex /c-plane. This support depends on the space of functions 
to which belongs the initial datum q{x,t), and, for q in the Shwartz space, we 
recall in the appendix the main results of the spectral theory for the Zakharov- 



Shabat spectral problem (5.11). There it is shown that the distribution R has 
the following structure 

, . _ ^f -a^(fc,t)(5-(fc/) 

'^' '''" ^\a+{k,t)6+{ki) 

^N± ( Cn{t)^{k-k-)\ ^2ika:,x( 



+2.E„U(c+ (,),(,„,.) "^^0 -Me^^'^^-^n^.U) 



when we chose two basic solutions i^^{k, x, t) of ( ^.11 ) according to the following 
behaviors for k real 

P+ -e-^^'^^a-/p-\ + f 1 0\ 

1//3+ J^^' '^{e^'''^a+ l) ^^■^^> 

1/13- \ _ (\ e-^'^^-a-X ,_„. 

It should be clear that the behaviors at a; = — oo are actually a consequence of 
the choice of the behaviors at a; = +oo. Here above the new coefficients /3^ are 
explicitely defined from the datum of R in the appendix, we simply need here to 
remark that, from the property ( ^.12 ), we have the so called unitarity relation 



a+a- +P+f3- = 1. (5.17) 

The main constraint resulting from the above structure of the distribution 



i?, imposed on M to ensure self consistency of the equation (3.4), read 

m±(fc, t) = TO±(fc, t)5^{ki) + ^ m^{t)5{k - fc±). (5.18) 

No such constraint applies on 17 because it is a fuction multiplying R. 
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ii) Natural constraints 



Since we are interested on evolutions (5.1) where the integral on the r.h.s runs 
on the real axis only, it is readily seen on (|5.13|) that it is necessary to set first 



m±(i) = 0. 



(5.19) 



Then it will be necessary that p — duj be proportional to the distributions 
(5^(fc/), hence 



p(k, t) == p^ikR, t)S+{ki) + po {kR, t)S [ki 



(5.20) 



Inserting now the set of constraints ( ^.18 )( 5yT9| )( ^.20 ) into the evolution 
( ^.13| ), we obtain the following system 



qt 



2i 



2i 



,+ /',,+ \2J2ikx 



2Po A*iiA*i2 + 2po Aiii-"i2 + "^0 (Mil) e ' '"^ - ^0(^*12) e 



(5.21) 



n = — 2p+^+ ^+ + 2Po ^^21^^22 + "^0 (M21) e ' "^ - ^0(^^22 

TT 



+ (,,+ \2^2ikx 



where of course the matrix /i = {/i^ } solves ( p. 11 ) , and where the notation /ij ■ 
is defined in (2/7). I t tur ns out (see appendix) that indeed the functions ju,^ do 
obey the behaviors ( 5.15| ) ( 5.16 ). 



iii) Reduction constraints 
It is convenient for physical applications to look for reduction of the above 6-field 
system to a simpler one (actually a 3-field system). This is done by assuming 
the following reduction 



for which from (5.11) 



r(x, t) = <yq{x, t), a = ± 



Mil ^ M22' M12 ^ '^M2ii 



and consequently using (5.14) in (3.1) 



oc- 



(5.22) 
(5.23) 

(5.24) 



The evolutions ( |3.4D together with (5.21) then imply the following constraints 
ml(k,t) = cnfig{k,t), Po(k,t) = -p^{k,t). (5.25) 

It is convenient to define lighter notations through 

Too = TO^, po = Pq, ol = ol^ , /3 = /3+, C„ = C+, fc„ = fc+. (5.26) 
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Finally, considering all the simplifications resulting from the above three 
stages of constraints, the system (5.21)(^|Tl) reads 



qt 



2i f 

■^ J-oo 



+ \2 2i\x 



(5.27) 



IJ'ii,x = qiJ'21, fJ'2i,x-2ikfi2i= crq^iu, fi ^2 = '^1^-21 ■ 

Since the relevant datum for the evolution of q{,t) is the function po{k,t), we 
need to compute uj{k,t) out of (5.20) 



u}{k,t) 



+ 00 






+00 



d\ 



X-k-iO 



Po{X,t), (5.28) 



that is, by using (B^), (5.26) and (UM) 



1 P /-J \ 

u){k,t) = -i(po+Po) PJ T r(Po -Po)- 

TT I \ — k 



(5.29) 



It results that uj is purely imaginary which makes the constraints (5.24) com- 
patible with the evolution (3.4). 



5.4 Evolution of the spectral transform 



To acheive the proof of integrability of the system (5.1) with arbitrary boundary 
values for the fields (01,02), it is necessary to find the transformation from 
( p^ ) to ( 5.27 ) by expanding the vector (01,02) on a suitable basis constructed 
from the vectors (/in, /i2i)- To do that it is enough to compare their asymtotic 
boundary values (they both solve the same first order differential linear syst em) , 
and hence, from ( ^.15| ) (5.16) we can deduce for each of the cases (5.2) (5.3) 

(5.30) 

(5.31) 

(5.32) 

(5.33) 



0-1 \ T - , T + 2ifca: 

= hfJ-i + ^2^26 

02 

02 j "/3+''i +/3-^2'' 



+ ^2ikx 



- ^^ ,,+ -1- T ,, + p2ifca: 



Next it is necessary to compute the product 01O2 appearing in the equation 
(^.l|) and to express it in terms of the only four products appearing in (5.27), 
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namely A'nA*i2i MriA'r2' [fJ-iiY e^'^'''^^ and (^^g)^^^*'''^. This is possible by re- 
peatedly using the Riemann-Hilbert problem for the function /i, as shown in 
the appendix. Actually we discover there that the real part of po can be set to 



zero. Once this is done, it is the sufficient to compare the equations (5.1) and 
( ^.27 ) and to equate the different coefficients of the above mentionned terms. 
The result reads: 



i) In the case (|5J) with (^.30|) 



po = -i-g{a\h\' + \l2V) 
mo 



yff[a/i/2-ia(a|/ip + |/2p)]. 



(5.34) 



ii) In the case (5^) with ( 5.31 ): 



po = -t-g[a\Ji\^{l + a\a\^) + \J2\^\P\^ + JiJ2a^ + JiJ2a(3], 



(5.35) 



77r - 1 1 - - 1 

mo = y.gMiJ2/3(l--a|«|2)--/3(a)2jiJ2 + -a|/3|2(a|Ji|2-|J2|2)]. 



in) In the case {pM with (5.32) 



.2. 1 + cr|Q;|- 



l-CT|a|' 



a 



a 



Po = -t-g[ia\Ki\'' + \K2\^)^-^^^+KiK2^ + KiK2^] 

ITT — — 1 

mo = —g[cykiK2(3{P)-^ + -a{a\K^\'' + \K2\'')] 



(5.36) 



iv) In the case (|5.5|) with (5.33) 



Po 



mo 






(5.37) 



^g[|(l - ^a\am,L2 - \jL1L2 + \a{a\L,\'' - IL^I^)]. 



With these values in hand, the evolution of the spectral transform defined 



in (5.14) is given by 



at = 2Lja — 2irao, 
kn,t = 0, 

Cn.t = 2Lj(fc„)C„, 



(5.38) 
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where uj{k,t) is given from po{k, t) in (^.29) 



It is convenient also, for future use, to obtain the evolutions of the coefficient 
(3 and of the quantity |a|^. From the definition of (3 given in the appendix, it 
is easy to derive 



ZITT 



dX 



i^wn 



X-k-iO l-cr|a|' 



and from the above evolution ( 5.38 ) 

{\af)t = 2i{fhoa - moa). 



(5.39) 



(5.40) 



Remember finally that, within the reduction ( 5.22| ), the unitarity relation (5.17) 
becomes 

a|ap + |/3|2 = 1. (5.41) 



5.5 Comments 



We have shown here above that the time evolution (3^) of the spectral transform 
can be computed explicitely in terms of the asymptotic boundary values ( p.2| )- 
( p^ ) in the two cases a — ±, when q(x, t) obeys the coupled system (5^). As can 
be seen in the preceding subsection, these evolutions can be quite complicated 
and possibly not explicitely solvable. In such a case the property of integrability 
of the system (pj) is weakened by the constraints of the boundary values. Still 
the spectral transform remains a quite usefull tool as it maps the problem of 



the solution of (5.1) (3 variables with one on the infinite line) into the simpler 
problem of solution of ( ^.38 ) (2 variables both on compact support). 

We have displayed in [pi a few applications of the present results for physical 
problems of resonant wave interaction. The important result to emphasize here 
is the strong dependance of the nature of the solution on the chosen boundary 
values. For instance, in the case of self-induced transparency, we obtain that 
|ap— >Oasi— >c» and this is precisely the transparency property. For stim- 
ulated Brillouin scattering in plasmas we get jap— >last— s-oo, and this is 
precisely the property of total reflexivity. We also give an example of special 
boundary values in stimulated Raman scattering for which jap — > oo as i ^ tg 
(finite), in that case the general solution blows up in finite time. There are other 
physi cal situations (under study now ||l^) for which the time evolutions ( 5.38 ) 
( ^.39| ) and ( 5.4C ) are not explicitely solvable. In that case, although the sys- 
tem possesses a Lax pair and a spectral transform, its integrability is partially 
destroyed by the constraints due to the boundary values. For instance, such a 



situation would be encountered here if considering the evolutions (5.38) (5.39) 
and (5.4C) in the case ( 5.35| ) with Ji J2 ^ and Ji external (not functionals of 
a,/3...). 
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As a last comment it is instructive to consider the case of self-induced trans- 
parency which corresponds to the choice |^ 

(T = -, R'l = 1, K2 = (5.42) 



in (|5.36| ). Hence we get 



TT 1 — lap n 

Po = '^-^9 , , I 12 ^ mo = -i-ga, 5.43) 

8 1-1- \ar 4 



and the evolution (5.3S) reads 



n 1 „/■ dXg l-|a|2 

-ga--aP - — ^ '—{-. (5.44) 



2 / A-fcl 



We have also from (5.40) 



{\a\')t^-ng\a\', (5.45) 

which proves the transparency: |ap vanishes exponentially as t ^ 00. The 
interesting (curious) fact is that, originally in ||ll|, the time evolution of the 
spectral transform has been written for the quantity 

p^aaP{/3)-\ (5.46) 

which is, in the language of the scattering theory, the r eflect ion co efficie nt to 
the left. It turns out that, by using successive^ 
evolution of p{k, t) can be written with help of 



the left. It turns out that, by using successivel y (|5.44| ) ( ^.45 ) and ( 5.39 ), the 



"-4/13^ "«■ <"" 



The point is that this evolution is Zmear while the equivalent one for a (5.44) 
is not. This property lead the people thinking that the situation for coupled 
systems was much the same as that of scalar systems, namely that the difference 
concerns only the dispersion relation which from polynomial becomes a singular 
function of k. This has produced a lot of works § [0 Q [|l|(l|l [|l8| , including 
those of the author Q [g| Q , where the linearity of the evolution of the spectral 
transform has been taken a priori, which implicitely imposes constraints on the 
boundary values. We see here that the linearity of the evolution can occur in 
particular situations, but is never a general property. The first instance of a 
system of coupled waves corresponding to a nonlinear evolution of the spectral 
transform has been given by KAUP in [^ in the case when the spectral problem 
is the Schrodinger scattering problem. This problem has then been further 
studied in pOJ and an extension to arbitrary boundary values is also feasible 
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6 Appendix: spectral and d problems 

6.1 Spectral problem and basic solutions 

We briefly recall here the basic results on the Zakharov-Shabat spectral problem 
which we write for the 2x2 matrix /i(fc, x, t) 

Its two fundamental solutions, say fi , are determined by (a t-dependence is 
understood everywhere) 



1, ^x^^{k, x)) \l)^\ -C d^r{Oti^^{k, 



(A.2) 
(A.3) 

(A.4) 
(A.5) 



The first column vector nf of the matrix /i+ is meromorphic in Im (fc) > 
where it has a finite number N'^ of poles k^ (assumed to be simple). The 
second vector fi2 is holomorphic in Im (fc) > 0. The vector ^^ is holomorphic 
in Im (fc) < while the second one y^^ is meromorphic in Im (fc) < where 
it has a finite number N~ of pole k~ (simple). One can check directly on the 
above integral equations that we have the relations 

Resn+{k)^tC+fi+{k+), (A.6) 



Res ^i^ik) = -iC^nlik-), (A.7) 

which define the normalization coefficients C„ . 

The function /i(fc) defined as /i+ in the upper half plane and fi^ in the lower 

is then discontinuous on the real fc-axis. Its discontinuity can be expressed 
simply in terms of fi itself as 

M+ - A^r = e'''=""+M^, (A.8) 

fi+ -fi^ ^-e-^'''^a-fi^, (A.9) 
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which define the reflection coefficients a^{k): 

r'+oo 



a+(fc) 



d^mt^tiiKOe 



-2ike, 



+ CXD 



{k)= dU{0f^22ik,0e'''^- 



We can also define the transmission coefficients P^ 

f-\-oo 



/ + OC 
-oo 



(A.IO) 
(A.ll) 

(A.12) 
(A.13) 



The integral equations (|A.2|) to (A. 5) give the foUowing behaviors at large x: 



13+ -e-2»fe^a-//3- 
1//3+ 



1 







-oo V e2»'=^a+ 1 



l/[i- Q \ ^ ^_ J I e-2«'=^a- 



^2^kx^+|P+ p- J _^>^ +^ I 



(A.14) 



(A.15) 



where we have used that from (A.l), det(/i) is z-independent. Consequently we 
have also the unitarity relation 



a+a- +P+(3- = 1 



(A.16) 



To obtain the asymptotic behaviors — e~^'^^a~//3 ~ oi nt-? ande~^**''^a+//3+ of 
IMjas X go es to — cx), one must use the relations ( A.8 ) and (|A.9| ) repeatedly in 
(A. 2) and ( A. 4 ) respectively. 

Solving the direct scattering problem consists in solving for given Q{x) the 
integral equations (A. 2) to ([A.q ) and then calculating the spectral data S: 

S = {a±(fc), A; e R; fc±, Ct,n = 1, .., 7V±}. (A.17) 



6.2 (9-problem and spectral transform 

Since the work |p2|, we know that the solution of the inverse problem, i.e. the 
reconstruction of Q from S, is given by a Cauchy-Green integral equation which 
solves a 9-problem for the function /i(fc) previously defined. Actually the d- 
problem is simply the formula which summarizes the analytical properties of 
/x(fc), it reads 

^^(fc) = M(fc)^(fc) (A.18) 
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where the distribution R{k) is the spectral transform and is given from the 
spectral data S by 

u(j,\-^f -a-ik)S-{ki)\2tkasx 

^^'''~ 2\a+{k)S+{ki) le - 



N^ 



+ 27r V ( _^ ,^ ^, Q ^(fc - ^n ) ] ^2tka3X 



C+S{k - k+) 







This formula can be demonstrated simply by noting that 
fceR: d^i{k)= ^fi{k)[d+(kj)-d-{kj)] 



(A.19) 



Im (fc) > : dn{k) 



-2iTTY,6{k - k+) Res fi+(k) , 



(A.20) 



lm(fc)<0: dn{k)^ ( , 2i7r2(5(fc - fc-) Res^^(fc) 

k„ 



The above 9-problem is completed by the asymptotic behavior of fi at large k, 
which is obtained from (A. 2) to (A. 5) by integration by parts and reads 



liik) = l + 0[l 



(A.21) 



Finally the solution of the equation (A. 18) obeying the above behavior is ob- 
tained by solving the following Cauchy-Grccn integral equation: 



,l^ - 1 ffdXAdX ,,,„,,, 



(A.22) 



By comparaison of the different powers of 1/fc in the equation (A.l), we 
readily obtain 

Q^i[c73,^i^% (A.23) 

where we have defined /i'^' through 



i{k) = 1 + ^2^' 



,(") 



(A.24) 



Therefore the inverse problem is solved by the integral equation (A.22) 
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